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Abstract 

We investigate the formation of trapped surfaces in cosmological spacetimes, 
using constant mean curvature slicing. Quantitative criteria for the forma- 
tion of trapped surfaces demonstrate that cosmological regions enclosed by 
trapped surfaces may have matter density exceeding significantly the back- 
ground matter density of the flat and homogeneous cosmological model. 
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I. INTRODUCTION 



In our previous work ( [[I]]) we investigated the formation of trapped surfaces in various 
cosmological models. Recently we have found a particularly useful formulation ( 0) of 
the spherically symmetric Einstein constraint equations that allowed us to improve our 
early estimates ( ||) for conditions determining the appearance of trapped surfaces. In 
the present paper we apply the new formalism to spherically symmetric cosmologies. As a 
result we find stronger criteria in spacetimes that were investigated previously and, more 
importantly, we are able to deal with hyperbolic universes where our previous attempts have 
failed. 

The order of the article is as follows. The first section presents the formalism. In Section 
2 we deal with the main results. Section 3 shows that regions enclosed by trapped surfaces 
must be invisible to external observers. The last Section contains conclusions of which 
the most important is that energy density inside cosmological regions enclosed by trapped 
surfaces may exceed significantly the average energy density of the flat and homogeneous 
cosmological model. 

There exist three homogeneous spherically symmetric cosmologies. The three are 

i) the closed (k=l) cosmology with metric 

ds 2 = -dr 2 + a 2 (r)[dr 2 + sin 2 rdtt 2 }, (1) 

ii) the open flat (k=0) cosmology with metric 

ds 2 = -dr 2 + a 2 (r) [dr 2 + r 2 dQ 2 ] , (2) 

iii) 

ds 2 = - dr 2 + a 2 (r) [dr 2 + sinh 2 rdtt 2 } , (3) 

where dfl 2 = dO 2 + sin 2 9d(j) 2 is the standard line element on the unit sphere with with the 
angle variables < (ft < 2n and < 9 < n. 
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The geometric part of the initial data set of the Einstein equations consists of the intrinsic 
three-geometry and the extrinsic curvature K ab which is essentially the first time derivative 
of the metric, all given at some time (say r = 0). The intrinsic geometries are respectively 

a 2 (r)[dr 2 + sin 2 rd£l% (4) 

a 2 (r)[dr 2 + r 2 dQ 2 }, (5) 

a 2 {T)[dr 2 + smh 2 rdtt 2 }, (6) 

and in each case the extrinsic curvature is pure trace 

K ab = Hg ab (7) 

where if is a time dependent function that is constant on each slice r = const. It is called 
the Hubble constant and it is given by H = 

In the general case initial data consist of the quartet (gfy, K^, p, J$) where gij is the 
intrinsic metric, Kij is the extrinsic curvature, p is the matter energy density and J; is the 
matter current density. These cannot be given arbitrarily but must satisfy the constraints 

{3) TZ - KijK ij + (trK) 2 = 16vrp (8) 

ViK lj - VHrK = -8nJ j (9) 

where ^7Z is the scalar curvature of the intrinsic metric. 

The momentum constraint, (§), is identically satisfied in the case of homogeneous cos- 
mologies (with Jj = 0) and the hamiltonian constraint, @, reduces to 

lQ 7T p = ^l + QH 2 (10) 
a 2 

where k is 1, 0, —1 in the closed, flat and hyperbolic cosmologies, respectively. Thus we can 
conclude that all slices of the constant coordinate time have a uniform energy density p 
which is at rest. 



In this article we wish to consider data for spherically symmetric cosmologies which 
either in the large approximate the standard cosmologies or asymptotically approach them. 
In all cases we will make the assumption that the initial slice is chosen so that the trace of 
the extrinsic curvature is constant on the slice. In order to retain the link with homogeneous 
cosmologies we define trK = 3H. 

The initial data we consider is a spherically symetric set consisting of a three-metric 

dsf 3) = a 2 dr 2 + b 2 (r)f(r)dn 2 , (11) 

an extrinsic curvature 

K r r =H + K(r), K$ = H-K/2, = H - K/2 (12) 

an energy density p{r) and a current density The function /(r) will be one of the set 
sin(r),r, sinh(r), depending on the type of cosmology. 

There are some useful geometric quantities that can be defined. One of them is the 
proper distance from the center of symmetry given by dl = adr. The Schwarzschild (areal) 
radius R is given by R = bf. The mean curvature of a centered two-sphere as embedded in 
an initial three dimensional hypersuface is 

2dlR n*\ 

P= nr - (13) 

In a general spacetime we may investigate the geometry by considering the propagation 
of various beams of lightrays through a space-time. These beams in general will shear and 
either expand or contract; a number of (optical) functions will be required to describe their 
propagation. In a spherically symmetric spacetime we focus our attention to light rays 
moving orthogonally to two-spheres centered around a center of symmetry. We need only 
two functions. These are the divergence of future directed light rays 

e = h ut^- R (14) 

ti CLT out 

and the divergence of past directed light rays 



where -^ out is the derivative along future-pointing outgoing radial null rays and 4p. is the 
derivative along future-pointing ingoing radial null rays. One interesting property of 9 and 
9' is that they can be expressed purely in terms of initial data on a spacelike slice. In the 
spherically symmetric case we have 

9 = p - K r r +trK = p- K + 2H, (16) 

and 

9' =p + K; -trK =p + K -2H. (17) 

This means that 9 and 9' are three-dimensional scalars. They are not four-scalars, since 
they depend on a choice of affine parameters along the null rays. 

In the homogeneous universes we find that pR = 2, pR = 2 cos(r) and pR = 2cosh{r) in 
the k — 0, 1, — 1 cases respectively and 



R9 = 2 + 2RH = 2 + 2\ i^-R, R9' = 2 - 2RH (18) 



for k=0, 



R9 = 2 cos(r) + 2RH = 2 cos(r) + 2aH sin(r) = 2 cos(r) + 2\ ( ^ - 1) sin(r) 



R9' = 2 cos(r) - 2\ (^Ml _ i) sin ( r ) (19) 



for k=l, 



R9 = 2 cosh(r) + 2RH = 2 cosh(r) + 2\ ( np ° a + 1) sinh(r), 



R9' = 2 cosh(r) - 2^ (^Ml + i) sinh (r) (20) 

for k=-l. 

A surface on which 9 is negative is called, after Penrose ||), a future trapped surface 
and a surface on which 9' is negative is called a past trapped surface. The occurrence of 



such surfaces in a spacetime is an indication of the fact that the gravitational collapse is 
well advanced. In the case of homogeneous closed cosmologies future trapped surfaces exist 
for any r > cot -1 (aH). In neither k=0 nor k=-l is R9 ever negative if H > 0. 

In this article we consider a universe that is homogeneous in the large but that it is 
dotted with numerous spherical inhomogeneities, far from each the metric approaches the 
background metric of a homogeneous universe. If we center our coordinate system at a 
particular lump we expect that optical scalars approach the values given in (|18| , [19], |2(J) far 
away from the lump. In the case of closed cosmologies this limiting value is expected to be 
met for values of the coordinate radius r much less than 7r/2. 

We assume local flatness at the origin, i. e., lim^Q R9 = \im.R_^ R9' = 2 al- 
though this condition can be relaxed to allow for a conical singularity there, i. e., 
< lim^o R6, lim^o Bff < 2. 

II. MAIN CALCULATIONS. 

The spherical initial data must satisfy the constraints, which read, in terms of functions 
9 and 9' 

d^OR) = -%irR{p-j) - -^-[2(9Rf -6MR-A- 129RHR] (21) 

AR 

di{&R) = -8nR(p + j) - -^-{2(6' Rf - ORO'R - 4 + 129RHR] (22) 

AR 

where j = ji is the radial component of the matter current density normalized so that 
J 2 = j jk- We can manipulate equations ( f2~l~D and (|22|) to obtain 

di{&R9R) = -87i(p(9'R + 9R)+j(9R-9'R)) - ^-[(9R9'R - A)(9'R + 9R)}. (23) 
^ ' 2R 

Let us now assume that the total matter satisfy the dominant energy condition, i. e., p > 
Assume that 9R9'R > 4 at a particular point. Consider first the situtation where both 
9R and 9'R are positive. Then (9'R + 9R) > (-9'R + 9R) and p{& R + 9R) +j(9' R + 9R) > 0. 
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This means that both terms of ( p3[) are nonpositive and the derivative of the product 8R6'R 
is negative. On the other hand, when both OR and 6'R are negative and their product is 
greater than 4, then p(9'R + 9R) +j(9'R + 9R) < and the first term in (p3|) is positive. The 
second term becomes also positive, so that di(9R9'R) > 0. Thus in both cases if 9R9'R > 4 
then di(9R9'R) ^ 0. 

Let us now consider the expressions for the product of the two scalars 9R9'R in each of 
the three homogeneous cosmologies. We get 

R9R9' = 4 - AR 2 H 2 (24) 

for k=0, 

R9R9' = 4 cos 2 (r) - 4R 2 H 2 (25) 

for k=l and 

R9R9' = 4 cosh 2 (r) -4(^^ + 1) sinh 2 (r) = 4([1 - ^^]sinh 2 (r)) (26) 

3 3 

for k=— 1. 

In each of these cases we have R9R9' = 4 at the origin and never more than 4. We are 
considering initial geometries that locally are flat and asymptotically approach the homoge- 
neous cosmologies, so that both at the origin and far from the center the product R9R9' does 
not exceed 4. If it were to achieve a maximal value greater than 4 somewhere in between, 
then its derivative would have to vanish; but that is excluded in the preceding analysis. 
Therefore we have proven: 

Lemma 1. Assume that matter satisfies the dominant energy condition and that spher- 
ical cosmological data are locally flat at the center and are asymptotic to any of standard 
homogeneous cosmological models. Then 

R9R9' < 4. 

Remarks: 
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i) The above statement is true for any regular slice, with arbitrary (i. e., nonconstant 
on a part of a slice) tri^, assuming that the slice is asymptotic to a homogeneous constant 
mean curvature slice. 

ii) It implies the positivity of the Hawking mass on a sphere centered around a symmetry 
center; 2M# = R(l — RdRd ' ) cannot become negative on a fixed slice. 

Lemma 1 holds true for all three cosmological models. 

The main issue that we will address in this paper is the question of the formation of 
trapped surfaces due to concentration of matter. The result will be obtained through a 
careful analysis of (|2i|) . What we do is multiply (^Tj) by R, use (|H|) and write the resulting 



equation in the following form 

di(6R 2 ) = -8vri? 2 (p - j) + 1 + ]-9R9'R - ^(9R) 2 + 39RHR. (27) 
The substitution of fllTf ) and (|I7|) into ( p7|) gives 

ott2 i 

d^R 2 ) = -8vri? 2 (p - j) + 1 + ~(pR + KR) 2 - R 2 K 2 + 2R 2 Hp (28) 

8tt 4 



or 



3k 1 

di{9R 2 ) = -8nR 2 (p - p + -j) + 2-{l--(pR+ KR) 2 ) - R 2 K 2 + ARHdiR, (29) 



where we used the relation ([10|) to eliminate the if 2 term and use the definition of mean 
curvature p. 

Let us integrate (p9|) from the origin out to a surface S. We identify 

r L ( s ) r 
AM = 4n R 2 (p-p )dl= dV(p-po) (30) 

as the excess matter inside a volume V(5') bounded by S and 

/•L(S) /• 

P = 4vr / i? 2 jc/Z = / d^' (31) 

JO JV(S) 

as the total radial momentum inside S. In this notation, the aforementioned integration 
yields 
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-2(AM - P) - ^- 2 V + 2L + ^ - dV(l - \(pR + KRf + R^) (32) 

where A is the area of the surface S and L is the geodesic distance of S from the centre. 
Below we will prove, in a series of lemmas, that under some conditions we can control the 
sign of the last integral. 

Lemma 2. Assume k=0, 1 cosmologies which are locally flat. If the energy condition 
P — Po — \j\ > — if^2 is satisfied out to an asymptotic region then 

2>\pR + KR\, 2>\pR-KR\. 

Lemma 3. In a data set that approaches the k = — 1 locally flat cosmology, if the energy 
condition p — p — \ j ' \ > is satisfied inside a sphere S then 

2>(pR + KR), 2>{pR-KR). 

Before proving the two lemmas, let us formulate two main results that give sufficient 
conditions for the formation of trapped surfaces. 

Theorem 1. Given data which approaches either the k=0 or the k=l locally flat cos- 
mology, if the energy condition p — p — \ j : \ > is satisfied out to an asymptotic region 
and if 



AM-P>--^y + L+^ (33) 



at a surface S then S is future trapped. 

Proof of theorem 1: the result follows directly from eq. (^) and the estimate of 
Lemma 2. 

Theorem 2. Assume that normally ingoing light light rays are everywhere convergent 
inside a volume V bounded by a surface S, 9' > 0. Given data which approaches the k = — 1 
locally flat cosmology, if the energy condition p — po — \j\ > is satisfied inside the volume 
V and if 
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Q IT A 

AM-P>-^-V + L + tp (34) 
87ra 2 47r 

at the surface S then there exists a surface inside S that is future trapped. 

Proof of Theorem 2. Assume that there is no future trapped surface inside S, i. e., 
9 = p + 2H — K > 0. Since we also assume that there is no past trapped surface, we may 
conclude that inside S p — K > —2H,p + K > 2H. We know that p is positive inside 
S because we have that p = (6 + 0')/2 and each of 9 and 6' is positive. We also have 
2 > pR — KR > —2HR and 2 > pR + KR > 2HR; the last inequalities follow from lemma 
3. If H > we have that pR + KR is positive and thus (p + K) 2 R 2 < 4 and the last integral 



of (^) is strictly negative. On the other hand, if if < we must have that pR — KR is 
positive and (pR — KR) 2 < 4 but we could have that pR + KR be negative. This can 
only happen while K is negative since we know that p is positive. In this case we write the 
integrand of (|3^ as 1 — \ {pR — KR) 2 —pKR 2 + K 2 R 2 . This is clearly nonnegative. Thus we 
also have in this case that the last term in ([32]) is negative. This contradicts the assumption 
that there is no trapped surface. Hence, under the assumptions of Theorem 2, there must 
exist a trapped surface inside S. 

In order to prove lemmas 2 and 3 we shall return to equations (f2T|) and (^3) and write 
them in terms of Rp, RK and RH. ( pl|) can written as 

d t (pR - KR) = -8nR(p -j-^-)-^(R p - RK) 2 + ±(Rp - RK)(Rp + RK) + i 

(35) 

and (122T) as 



d t (pR + KR) = -8nR(p + j-^-)- ^(Rp + RK) 2 + -^(Rp - RK)(Rp + RK) + 1. 

(36) 

We will prove first the upper bound of pR + KR, pR — KR, simultaneously for both Lemma 
2 and 3; this part of the proof does not depend on the type of a cosmological spacetime. 
Also, as it will become clear, the energy condition shall be imposed only inside a sphere S if 
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we are interested in finding the estimate inside S (as opposed to the estimations from below 
that require the global assumption made in Lemma 2). According to the conditions made 



in lemmas, the first term of either equation (35) or (Rq) is nonpositive. We show that in the 



situation of interest the remainders of each of the equations are also nonpositive. 

At the center of symmetry the quantities pR+KR, pR—KR are equal to 2, for all types of 
cosmology. This means that right hand sides of either (^) or (J36Q must be nonpositive and 
that the quantities in question start from the origin with the value 2 and start to decrease 
as soon as they meet either positive p + j — or p — j — ^— . 

Let us assume that further out one of the two, say pR + KR, rises up to 2 with pR — KR 
lagging behind. In this case we can write the non-material part of the right hand side of 
(j36|) as follows 

- ~(Rp + RK) 2 + ^{Rp - RK)(Rp + RK) + l = -l + -(Fp- RK) < 0. (37) 



Because the material part of ( p6|) is nonpositive, we get that di(Rp + RK) < so that 
pR + KR cannot exceed 2. A similar argument can be made for pR — KR. Thus Lemma 3 
and the upper bound of Lemma 2 are proven; as is clear from the above derivation, in order 
to have a bound that is valid inside a sphere S we need the energy condition that is imposed 
only inside S. 

The same reasoning can be applied to complete the proof of Lemma 2. We will show, 
that if one of the two quantities in question reaches the value -2, then at least one of them 
must be less than -2, thus breaking either the demand of geometries being asymptotic to a 
homogeneous cosmology in the sense expressed in equations (|18D and (|T9|). 

In order to show this we need the global energy condition of Lemma 2. Let us assume 
that there exists a point where pR + KR = —2, with pR — KR > pR + KR. Then the 



nonmaterial part of Eq. ( p6|) reads 



- -[2(Rp + RK) 2 + ^(Rp- RK)(Rp + RK) + 1 < (38) 

Eq. ( |3T)D implies now (assuming the energy condition) that pR + KR has to become more 
negative, if pR + KR < pR — KR and may stay at -2 in the case of equality only if the matter 
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contribution exactly cancels. However, if we can impose an outer boundary condition such 
that pR + KR > —2 then we get a contradiction. A similar argument works for pR — KR. 
The outer boundary condition is guaranteed in the cases of interest. Cosmo logical spacetime 
dotted with inhomogeneities have the property that asymptotically pR + KR and pR — KR 
approach values given by flT8p and ( |T9"D which must be strictly bigger than -2. That ends 
the proof of Lemma 2. 

It is interesting that we obtain an exact criterion with the constant 1; this suggests 
that the above theorem constitute a part of a more complex true statement that can be 
formulated for general nonspherical spacetimes. It suggests also that M(S) is a sensible 
measure of the energy of a gravitational system that might appear as a part of a quasilocal 
energy measure in nonspherical systems. 

It is clear that the analysis performed here can include cases where the sources are 
distributions rather than classical functions; in particular, we have no difficulty with shells 
of matter. All we get on crossing the shell is a downward step in 9 and 9'. More interestingly, 
we can extend the analysis to include conical singularities at the origin ( |]]), in a way 
analogous to that described in ( M). 



III. CONFINING PROPERTY OF TRAPPED SURFACES. 

In this Section we show that a region enclosed by trapped surfaces cannot be seen by 
external observers. This fact has been proven (without referring to the Cosmic Censor 
Hypothesis) by Israel ( H) . Here we will present a different version of the proof that is 
based on a 1+3 decomposition of a spacetime (as opposed to the proof of Israel, who used 
a 2+2 decomposition). 

We need the evolution part of the Einstein equations and the lapse equation. These are 

d t (5K r r - 2H(t)) = ^(5K r r ) 2 - ^ - -7=<9 r a + ^ + 8vr«T; + 3aH 2 - 3HSK r r (39) 
4 4 y a R 

and 
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V^a = a(^(^;) 2 4vr(p + 7?) + 3H 2 ) + 3d t H. (40) 

In addition we need the evolution equation of the mean curvature p of centred spheres 

d tP = ^(-5K r r + 2H) + 8na^= + ^(5K r r - 2H) (41) 

Using these equations we can find the full time derivative of 6 along a trajectory of null 
geodesies normal to centred spheres 

(d t + -^=)6 = ^6'6 - 8vra(-2^ + p + T r ) - a6 2 + 3aH6. (42) 
\/a \/a V<2 

Take now an apparent horizon, i. e., a centred sphere S of vanishing 6(S); fl4"2] ) implies that 
photons that start from S will forever remain inside an apparent horizon, if the strong energy 
condition — 2-^ + p + T r r > is assumed. Hence apparent horizons move faster than light 
in cosmological spacetimes (in contrast with asymptotically flat spacetimes, where they can 
eventually stabilize to the speed of light); they act as one-way membranes for non-tachyonic 
matter. This means that outside observers cannot detect any information from any inside 
region that is enclosed by a trapped surface. The only way to draw any conclusions about 
a piece of a spacetime that is enclosed by a trapped surface is through the observation of 
"long-wave" effects - through the attractive force that large massive objects exert on their 
surrounding. 



IV. DISCUSSION. 

Cosmological trapped surfaces that we discuss in preceding sections can, if they exist, 
accumulate an enormous amount of energy. Typically, as we have shown, the matter content 
of a trapped surface having a geodesic radius L is of the order L plus the background 
energy M# = 3I ^ n v (we neglect here the effects related to the possibility of a nonzero 
curvature of the space-like slice and the surface term Assume that there exists a 

trapped surface with a proper radius of the order of 1000 megaparsecs. Then its excess 
energy is of the order of 1000 (in units of megaparsecs). The present value of the Hubble 

13 



constant is about 50 — or (in units in which the speed of light c=l) 7^ — . 

s*megaparsec \ 1 ° > bwumegaparsec 

Therefore the expected value of the background energy inside the above ball is of the order 
0.5 * (■^^) 2 megapar sec = 14magaparsec, which is about 10 2 times less than the energy 
content that is needed in order to form, say a spherical massive shell that creates a trapped 
surface. We include this crude calculation just to point out that the formalism of general 
relativity does allow for cosmological regions with high concentrations of matter that are in 
principle invisible by external observers. 
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